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Self-Consistent Second Order Perturbation Theory 
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We apply self-consistent second order perturbation theory (SCSOPT) with respect to the on- 
site repulsive interaction U to study the Hubbard model in two dimensions. We investigate 
single particle properties of the model over the entire doping range at zero temperature. It is 
shown that as doping decreases toward half-filling cj-mass enhancement factor increases, while 
fc-mass enhancement factor decreases. The increase in u-mass enhancement factor is larger 
than the decrease in fc-mass enhancement factor, so that total-mass is larger than that in the 
non-interacting case. When particle number density per unit cell n is given by 0.64 < n < 
1.0 interaction enhances anisotropy of the Fermi surface, whereas at lower densities n < 0.64 
interaction suppresses anisotropy of it. Due to the decrease in fc-mass enhancement factor the 
density of states (DOS) at the Fermi level is suppressed. It is possible to understand the results 
within the framework of the weak coupling Fermi liquid theory. 

KEYWORDS: Hubbard model, two dimensions, Fermi liquid, self consistent, second order perturbation, Fermi 
surface, fc-mass, u;-mass 



§1. Introduction 

Since the discoAjery of the high-T c superconducting 
cuprates (HTSC), EP two dimensional layered materials 
have been attracting the attention of many researchers. 
Two dimensional Hubbard model is one of the promis- 
ing model to understand such materials. Investigations 
of the Hubbard model have a long, history and the un- 
derstanding of the one dimensionalo* and the infinite di- 
mensional cases tP is in a higher level. In the case of one 
dimension DOS at the Fermi level becomes zero with in- 
finitesimal interaction U. On the other hand in the case 
of infinite dimensions DOS at the Fermi level doesn't 
change at all until metal-insulator transition happens 
due to the locality of the self-energy. For the case of 
two dimensions the feature of the model hasn't been suf- 
ficiently elucidated. 

We use here the self-consistent second order perturba- 
tion theory (SCSOPT) with respect to the on-site repul- 
sive interaction U to study the Hubbard model in two 
dimensions. The SCSOPT satisfies the Luttinger sum 
rule, which is necessary to construct the Fermi liquid 
theoryJS The SCSOPT was formerly used for studying 
the infinite dimensional case© and was known to de- 
scribe the low energy excitations of the paramagnetic 
metallic phase of the model correctly in the weak cou- 
pling case. SB Therefore in the present paper we restrict 
our attention to the metallic phase of the model with- 
out magnetic order, which is suggested for a non-half- 



filling case in Quantum Monte Carlo simulations! 
On the contrary to infinite dimensional case the self- 
energy of the two dimensional Hubbard model has fc- 
dependence. The fc-dependence of the self-energy is con- 
sidered to play an important role in Mott transition in 
two dimensional systems such as Cai^xSrxVOaB' and 
HTSOEP The fc-dependence effects of the two dimen- 
sional Hubbard model were examij^cLsc.far by second 
order perturbation theory (SOPT) BB0 Jiiuithe fluc- 
tuation exchange (FLEX) approximation.ollS' In the 
case of SOPT it was reported that for sufficiently large 
U n(fi) exceeds unity in spite of /x < around \x ~ 0, 
which leads to unphysical negative charge susceptibil- 
ity.EJ Then we define U\ (jj) so that in the region U > U\ 
SOPT is unphysical. In addition SOPT satisfies the Lut- 
tinger sum rule only approximately for small U except for 
half-filling. This disagreement between the Fermi surface 
volume and the expectation value of the number operator 
is also remarkable near the half-filling. Then we define 
£/ 2 so that in the region U > U2 physical features are af- 
fected by that disagreement. We can use SOPT only in 
the region U < min{Ui, 1/2}- There are not such limita- 
tions in SCSOPT. Meanwhile the FLEX approximation 
includes SCSOPT diagrams, but it would be meaningful 
to distinguish the contribution from SCSOPT diagrams 
from the contribution from the other diagrams. 

We adopt the unit system in which h = 1 and fcs = 1. 
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§2. Model and Calculations 

The Hubbard Hamiltonian is 

H = ~ l E ( C *W + h.c.) + nu^ , (2.1) 

n.n.,a i 

where t is the transfer integral between the nearest neigh- 
bor sites and U stands for the on-site repulsive interac- 
tion. 

The one paxticle Green's function at zero temperature 
is defined asE^P 

G Ra {t) = -i < T^(fi,t)V4(O,0) > , (2.2) 

where ip a (R,t) denotes a particle field of electrons and 
(...) and T mean thermodynamic average over grand 
canonical ensemble and a time ordering operator, respec- 
tively. 

The Green's function in (fe, w) space can be expressed 
by the self-energy as 

( k =e k -» (2.3b) 



where 



and 



e k = —2cos{k x ) — 2cos(k y ) 



(2.3c) 



We define here the retarded self-energy and the retarded 
Green's function as 



and 



yWr.A - J E fcM UJ>0 
h k {UJ) ~ \ S* (w) u < 



(2.4) 



(2.5) 



The Green's function can be expressed by the spectral 
function as 



d£- 



where 



lj — E — i<5sign(i?) 
1. 



(2.6) 



Then it follows 



G ka {t)=\0{-t) f dEe- iEt+St A ka (E) 

J - — OC 

poo 

-iO{t) dEe- iEt - 5t A ka (E) (2. 
Jo 



-t)a k Jt)-i0(t)b ka (t) . (2.9) 



Here we define 



a ka (t) = f dEe- iEt + St A ka (E) (2.10a) 

J — OO 



and 



b ka it) S / dEe-^A ka {E) 
Jo 



(2.10b) 



Real space Green's functions can be calculated as 



k 



The contribution of the Hartree diagram (Fig. 1) to 
the self- energy is 

y-, (Hartree) , 



(2.12) 



As we consider the paramagnetic state in the present 
paper, 



£, (Hartree) ^ = ^ 



J ka 



(2.13) 



This static contribution to the self-energy can be incor- 
porated into the chemical potential, and as a result the 
redefined chemical potential takes the value at half- 
filling. 




Fig. 1. Hartree diagram. The wavy and double solid lines repre- 
sent the on-site repulsive interaction and the exact one particle 
Green's function, respectively. 



The second order contribution (Fig. 2) to the self- 
energy is 

E^(t) = U 2 G^(t)G^_(t)G^„_(-t) , (2.14) 



Ra^ ' R- 



-R 



where G R (t) is the Green's function of the non- 
interacting system. Similarly the SCSOPT contribution 
(Fig. 3) to the self-energy is 

(2.7) E<£(t) = U 2 G% c Ht)G% c Ut)G%_(-t) , (2.15) 



R(T 

where 



(2.16) 



In real space the retarded self-energy of SCSOPT is ex- 
pressed as 

Sgf c \t) = ^e{t){b Ra {t)b R _ a {t)a_ R _ a {-t) 
+ aR,(^R.At)b_ R _ a (-t)} 

= U 2 0(t){G R ^)G R c l(t)Gt c R _ a (-t) 
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G 



(sc)* 
-Ra 



(*)} 



(2.17) 



The real and the imaginary parts of the retarded self- 
energy satisfy the Kramers-Kronig (K.K.) relation, 



ReE 



Ka 7T J 



ImEi 
dw' h 



(2.18) 



The algorithm for numerical calculations is as follows. 

(1) We choose the initial E^ (ui), for example, the self- 
energy of SOPT. 

(2) We calculate G Ra (t) using eqs. @, ([Tj]), (2.8) 
and (2.11). At this time, since Au(u>) is localized 
with respect to energy within the extent of the band 
width, the integration in eq. (2.8) can be reduced to 
an appropriate finite energy interval. 

(3) Using eq. (2.17) we calculate En' (t) which goes to 
zero more rapidly than Ed (i) as t increases. 

(4) We calculate EV ,'(uj) using the Fast Fourier Trans- 
formation (FFT) . At this time it is efficient to calcu- 
late first ImEjk ( w ) which is localized and calculate 
ReE^(cj) using the K.K. relation eq. fl2.18| ). 

(5) If we iterate this procedure until convergence occurs, 
we obtain the SCSOPT sL R) (w). 

Here we present the results for a lattice with 256 x 256 
sites and the anti-periodic boundary condition. The 
FFT from time to energy variable was performed on 512 
points. The energy spacing 5uj is 0.06 and the momen- 
tum spacing Sk = 27r/256 is 0.025. All of the numerical 
results are calculated for zero temperature. A typical 
run to solve the Hubbard model with a fixed chemical 
potential value takes 1 hour on a DEC alpha station 250 
4/266. 

Hereafter we drop the spin indices for physical quan- 
tities as we are in the paramagnetic phase. 



£(R,t) 
G(R,t) 




= — > — + -►- 



Fig. 3. Diagrammatic representation of SCSOPT. The solid, the 
double solid, and the wavy lines represent free Green's func- 
tions, renormalized Green's functions, and on-site repulsive in- 
teractions, respectively. 



3. 1 Chemical potential shift 

As a function of density the interaction-induced chem- 
ical potential shifts are shown in Fig. 4. The chemical 
potential doesn't shift for half- filling (n = 1) case due to 
particle-hole symmetry, and has a maximal shift —0.261 
at n = 0.47 for U = 4. 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 



density 




Fig. 4. Chemical potential shifts v.s. particle density, fi and [iq 
represent the chemical potential of the interacting and the non- 
interacting system, respectively. 




Fig. 2. Diagrammatic representation of the self-energy in second 
order in U. The wavy and the solid lines represent on-site repul- 
sive interactions and free Green's functions, respectively. 



§3. Results 



In the weak coupling theory the chemical potential 
shift is related with the w =pJ3 self-energy on the non- 
interacting Fermi surface as EZP 



fyt~<ReE fe (0) >fso 



(3.1) 



where k p Q is the Fermi momentum of the non- interacting 
system and < ... >fso means the average over the non- 
interacting Fermi surface 

<r Rpv rm >> - ^k^k-fo^k (0) 

< Kelj k F S u ) >fso= 17 ^ ' 

(3.2) 



where /xo is the unperturbed chemical potential. 

Figure 5 shows the SCSOPT ReE feFo (0) as a func- 
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tion of n. The curve of chemical potential shift is put 
between the curve of ReE^, F0 (0) °f C 71 ") 0) direction and 
that of (ir, tt) direction. The variation of ReE^ (0) on 
the non-interacting Fermi surface is much smaller than 
ReE^, (0) itself. Therefore we can conclude that the 



deviation of S/i from eq. (3.1) is small compared with 5/i 
itself. 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.6 

density 



Fermi surface, and nu is a unit vector normal to the 
non-interacting Fermi surface and always is directed out- 
wards from kpo. <p represents the polar angles of posi- 
tions on the non-interacting Fermi surface in units of 
degree, which is described as Arctan(fc^oi//fcFOa;) mathe- 
matically. At densities 0.9, 0.8 and 0.7, which are larger 
than 0.64 Sk F is monotonously decrease as <j> changes 
from 0° to 45°. On the other hand at densities 0.5 and 
0.3, which are smaller than 0.64 Skp is monotonously in- 
crease in that region of tf>. Then there is a unique inter- 
section of interacting and non-interacting Fermi surfaces 
in that region of <fi except for the cases of n = 1 (due 
to particle-hole symmetry) and 0.64. That intersection 
moves from 0° to 22.5° monotonously as n changes from 
1 to 0. 

Then from Figs. 6 and 7 it is concluded that at den- 
sities 0.64 < n < 1.0 interaction enhances anisotropy 
of the Fermi surface, whereas at densities n < 0.64 it 
suppresses anisotropy. 

Fermi surface deformation is very small and the topol- 
ogy doesn't change. 



Fig. 5. £^ (0) of SCSOPT as a function of n for interaction 
strength U = 4. The results are shown for (tt, 0) and (it, it) 
directions on the non-interacting Fermi surface. Chemical po- 
tential shift for the same interaction strength is also shown for 
comparison. 



3.2 Fermi surface 

The excitation energy of the quasi-particle is deter- 
mined by the secular equation 



E 



k -Z k -ReX k (E k ) = 



(3.3) 



The Fermi surface of the interacting system is deter- 
mined by the condition E k = 0, hence 



£fc F +ReE feF (0) = 0, 



(3.4) 



where kp is the Fermi momentum of the interacting sys- 
tem. 

As shown in Fig. 6 at densities 0.64 < n < 1.0 the 
Fermi momentum of the interacting system is larger than 
that of the non-interacting system in (tt, 0) direction, 
whereas in (tt, tt) direction it is smaller than that of the 
non-interacting system. On the other hand at densities 
n < 0.64 the Fermi momentum of the interacting sys- 
tem is smaller than that of the non-interacting system 
in (tt, 0) direction, whereas in (tt, tt) direction it is larger 
than that of the non-interacting system. 

Figure 7 shows the Fermi surface deformation for var- 
ious densities. Skp is defined as 



Sk F (k F0 )=(k F - k F0 ) ■ n k 



(3.5) 



where we take the non- interacting Fermi momentum k F q 
so that k F — k F Q is perpendicular to the non- interacting 



(7C, o) 

direction 



(71,71) 

direction 



).0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 

density 



Fig. 6. Fermi momentum shifts induced by the interaction over 
the entire doping range in SCSOPT. kp and kpo denote the 
Fermi momentum of the interacting and the non-interacting sys- 
tem, respectively. The results are shown for (tt, 0) and (tt, it) 
directions on the Fermi surface. Interaction strength (7 = 4. 



In the weak coupling case eq. (3.4) is reduced tdL 
^-ReE fepo (0) 



Skp(k 



FQ) 



(3.6) 



where u?„ = |Ve^. | is the Fermi velocity of the non- 



interacting system. Putting eq. (3J) into eq. ( p.6[) we 

get 

< ReEi, (0) > FS0 -RcEi, (0) 
Sk F (kpo) ~ kF0 n =£aii , (3. 7 

Halboth and Metzner calculated the second order co- 
efficient of Sk F in U using eq. K]\) in SOPTO The 
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Fig. 7. Fermi surface deformation induced by the interaction for 
various densities. The horizontal axis represents the polar an- 
gles of positions on the non-interacting Fermi surface in units of 
degree. Interaction strength U=4. 



results of Fig. 7 are almost the same as the results 
of Halboth and Metzner except for the scale of the 



y-axis. This means that in SCSOPT eqs. (fU|) and ( pTq ) 
are in fact satisfied approximately and low energy part 
< £<, 2) (0, M oH) >fso -Sj? (0,/ioW) used in their 
work is almost the same as < £/ (0,/i(n)) >fso 
_£,(sc) ^o, jLt(ri)) except for the scale. The difference in 
the scale originates in the fact that in SOPT Skp in- 
creases in proportion to U 2 , whereas in SCSOPT it isn't 
in proportion to U 2 . 

As shown in Fig. 5 at densities 0.64 < n < 1.0 
Re££, Fo (0) m i 11 ! 0) direction is smaller than that of 
(7r,7r) direction, whereas at low densities n < 0.64 it 
is larger than that of (ir, tt) direction. These features 
combined with eq. (3.7) account for the results in Fig. 6. 



3.3 k -mass and lo -mass 

The denominator of the Green's function is expanded 
around E k as 



, , (9Re£i.(u; 
-£ fe -Re£ fe M^(l 



j=E 



k 



(3.8) 

Here we define w-mass and w-mass enhancement fac- 
tors as 



7a 



1 



dReZ k (u) 



1 



where m is the mass in the case of U = 



(3.9) 



(3.10) 



We define fc-mass and fc-mass enhancement factors as 
kv 

— = |grad fe (e fe +ReE fe (0))|| fe=feF (3.11) 



Ik 



(3.12) 



Dispersion of the quasi-particle energy is expressed in 
terms of mass enhancement factors as 
kp 



E k = ( k ~ k F)-n k 



(3.13) 

where we take kp so that k — kp is perpendicular to the 
Fermi surface and nu is a unit vector normal to the 
Fermi surface and always is directed outwards from k p . 

Figure 8 shows that as n increases toward half-filling 
cj-mass enhancement factor increases, while fc-mass en- 
hancement factor decreases (Exceptionally 7^(90°) in- 
creases slightly in the region 0.81 < n < 0.87). The 
amounts of the increase in cj-mass enhancement factor 
and the decrease in fc-mass enhancement factor are larger 
in (0, 7r) direction than those in {it, it) direction. 

From eq. ( 3. 13 ) total mass enhancement factor is de- 
termined by 7^7fe. The increase in w-mass enhancement 
factor is larger than the decrease in fc-mass enhance- 
ment factor, so that total mass is larger than that of 
the non-interacting system. As shown in Fig. 9 total 
mass enhancement factor in (0, it) direction has a peak 
at n = 0.97, and decreases for n > 0.97 due to the de- 
crease in fc-mass enhancement factor. 

The decrease in fc-mass in this weak coupling the- 
ory is not sufficient to explain the experimental data of 
Mott-Hubbard^ystems such as Cai-xSixVOs,^ HTSC 
La2- x Sr x CuOjl$ and an organic conductor BEDT-TTF 
salt a-(ET) 2 KHg(SCN) 4 . @ 

Recently Yoda and Yamadat^P investigated the ef- 
fect of the long-range Coulomb interaction on the fc- 
dependence of the self-energy using the perturbation the- 
ory up to the second-order terms. For the half-filling case 
they found that the first-order term made a large reduc- 
tion of the fc-mass. .— . 

On the other hand OhkawaHP suggested that in the 
vicinity of the Mott transition the formation of light 
quasi-particle could occur due to intersite quantum spin 
fluctuations using 1 /d expansion method within the on- 
site Coulomb interaction. 

The relation between these theories and real systems 
is not clarified enough now. 



3.4 Density of States 



From eq. (3.8) 



z 



(3.14) 



{u-E k ) -ImS fc (w) ' 

If we assume the system can be described by the Fermi 
liquid theory, the coherent part of the spectral function 
is 



(3.15) 
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1.8 




0.2 0.4 0.6 0.8 1 

density 



Fig. 8. fc-mass and w-mass enhancement factors v.s. density of 
particles. 7^(90°) denotes the w-mass enhancement factor at 
the Fermi surface in (0, 7r) direction, and 7^(45°) denotes the 
one in (tt,tt) direction. 7^(90°) and 7^(45°) are the fc-mass en- 
hancement factors defined in the same way. Interaction strength 
(7 = 4. 



1.8 i- 

s_ 1.7 - 
O 

O 1.6 - 
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0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

density 



Fig. 9. Total mass enhancement factor v.s. density of particles. 
7^7*: (90°) denotes the total mass enhancement factor at the 
Fermi surface in (0, 7r) direction, and 7^7fe(45°) denotes that 
in (it, it) direction. Interaction strength U = 4. 



in momentum space 

D(0) = ~f k D u=0 (Q) . (3.18) 

As shown in Fig. 10 the peak width of DOS shrinks 
as interaction increases, and the peak position of the 
DOS moves toward the Fermi level. These features re- 
semble those in the infinite dimensional case.l$ On the 
other hand DOS at the Fermi level is suppressed with 
the decrease in the fc-mass enhancement factor. This 
feature is not found in the infinite dimensional case be- 
cause the self-energy does not have fc-dependence, and 
consequently DOS at the Fermi level remains the same 
as that of the non-interacting case.Ef 

Figure 11 shows the comparison of the DOS of SC- 
SOPT with that of SOPT for n = 1. In the case of 
SOPT the weight near the Fermi level decreases and two 
peaks build up at high energies. These two peaks axe re- 
lated with the upper and lower Hubbard bandsOliSH 
On the other hand the SCSOPT DOS does not have such 
two peaks. The lack of high-energy resonance peaks in 
SCSOPT indicates that the SCSOPT is a weak coupling 
theory and is not correct for large U. 




-5 -4 -3 -2 -1 1 2 3 4 5 6 7 8 9 10 



energy 



Fig. 10. Renormalized density of states versus energy for n = 0.8 
of the SCSOPT for interaction strengths (7 = and 4. Inset 
shows the magnification near the Fermi level. 



Then the jump in the momentum distribution function 
across the Fermi surface is given by Z^,^ . The density of 
states of the interacting system is defined as 



(3.16) 



ko 



Putting eqs. (|3.13|) and ( p.l5|) into eq. fl3.1q ), we find 
that density of states at the Fermi level is related with 
the fc-mass enhancement factor 

D (°) - 2 ?i / ds ^ . ( 3 - 17 ) 



(2tt)' 



FS 



where 



|Ve^ | and J FS ds denotes the linear inte- 
gration along the Fermi surface. If the system is isotropic 



3.5 Charge susceptibility 

If the chemical potential is represented by fit, charge 
susceptibility Xc is defined as dn/d/if Since we define /i 
as nt — Un/2 in § ^, charge susceptibility is described as 

dn 1 , 
Xc = ^-jj . (3.19) 

Charge susceptibility for U — is described as 

d/X (27r)2 7 FS V% 

For U = 0, Xc diverges logarithmically toward half-filling 
due to the Van-Hove singularity. 



X 
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Fig. 11. Renormalized density of states versus energy for n = 1 of 
the SCSOPT and the SOPT (second order perturbation theory) 
for interaction strength [7 = 4. U = case is also shown for 
comparison. 



In the interacting case dn/du 
in the Fermi liquid theory asE3 



in eq. ) is described 



<9Re£ ( , sc) (0) 
1 



Putting eq. fl3.15| ) into eq. ( ]3.21] ) we obtain 



dn 1 

dfi = 2 (2tt) 2 



where we define 



9ReS fe (0) 



(3.21) 



(3.22) 



(3.23) 



As shown in Fig. 12 the charge susceptibility Xc is 
suppressed by interaction near half-filling. This suppres- 
sion is caused by the factor 7 M in eq. (3.22), and can be 
thought as the precursor of the Mott insulator charge 
gap. Similar features have been reported in other theo- 
ries which are in the framework of the Fermi liquid the- 
ory, such as the theory of infinite dimensional Hubbard 
modelLy and second order perturbation theory for Xc of 
the two dimensional Hubbard modelH 

On the other hand our result is inconsistent with , 
rukawa and Imada's Quantum Monte Carlo studj 
where Xc diverges as n goes to 1. 

§4. Conclusion 

We have applied the self-consistent second order per- 
turbation theory to study the two dimensional Hubbard 
model, which is a weak coupling theory and is expected 
to describe metallic phase of the model correctly. The 
SCSOPT satisfies the Luttinger sum rule exactly, which 
is necessary to construct the Fermi liquid theory. 

First we have calculated the Fermi surface deforma- 
tion, and the result shows that the interaction enhances 
the anisotropy of the Fermi surface near half-filling and 



Fig. 12. Charge susceptibility Xc as a function of the particle 
density. 



reduces the anisotropy of it at lower particle densities. 

The peak width of DOS shrinks as interaction in- 
creases, and the peak position of the DOS moves toward 
the Fermi level at non-half-filling densities. These fea- 
tures resemble those in the infinite dimensional case. On 
the other hand DOS at the Fermi level is suppressed with 
the decrease in /c-mass enhancement factor. The decrease 
in A;-mass enhancement factor is small and about 10 per- 
cent for U = 4 at half-filling where it is maximum. This 
feature can not be found in the infinite dimensional case. 

In the weak coupling theory the chemical potential 
shift induced by the interaction is related with the av- 
erage of the lu = self-energy over the non-interacting 
Fermi surface, and the Fermi surface deformation is re- 
lated with the variation of it on that Fermi surface. 
These relations are well satisfied in SCSOPT for the two 
dimensional Hubbard model. 
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